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We investigate a Hamiltonian lattice version of the two-dimensional Wess-Zumino model by Quantum Monte 
Carlo simulations. In order to study the pattern of supersymmetry breaking, we measure the ground state energy 
and the correlation length along a trajectory approaching the continuum limit. The algorithm is very effective in 
measuring the ground state energy, and adequate for the correlation length. 



1. INTRODUCTION 

Numerical simulations of lattice field theories 
are usually performed in the Lagrangian formula- 
tion. Nonetheless, we think there are very good 
reasons to develop numerical simulation tech- 
niques for the Hamiltonian approach |l| : powerful 
many-body techniques are available M, which al- 
low the direct computation of the vacuum wave 
function properties; fermions are implemented di- 
rectly and need not be integrated out; properties 
like the mass spectrum are more immediate. Fi- 
nally, universality checks between the Lagrangian 
and the Hamiltonian formalism are very welcome. 

2. THE MODEL 

We study the Hamiltonian lattice version of the 
two-dimensional Wess-Zumino model described 
in Refs. [|[^; we only wish to highlight here the 
main features of the formulation. 

In the Hamiltonian formalism, since H is con- 
served, it is possible to preserve exactly a 1- 
dimensional subalgebra of the original supersym- 
metry algebra, i.e., we can write H — Q^, where 
Q is a fermionic charge. This subalgebra is 
enough to guarantee some of the most impor- 
tant property of supersymmetry, including a non- 
negative spectrum, and pairing of fermionic and 
bosonic states of nonzero energy; spontaneous 
breaking of supersymmetry is equivalent to a 
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strictly positive ground-state energy Eq; the full 
supersymmetry algebra is recovered in the con- 
tinuum limit together with Lorentz invariance. 

In order to obtain a Hamiltonian free of fermion 
sign problems, and therefore amenable to Quan- 
tum Monte Carlo methods, we adopt free bound- 
ary conditions, with lattice size L = 2 (mod 4). 

The model is parametrized by a prepotential 
V{(j)), an arbitrary polynomial in the bosonic 
field. The two-dimensional Wess-Zumino model 
is superrenormalizable; fields do not renormalize, 
and only V{(j)) needs to be normal ordered. 

In strong coupling at leading order, the model 
reduces to independent copies of supersymmetric 
quantum mechanics, one for each site; supersym- 
metry is broken if and only if the degree of the 
prepotential V is even In weak coupling, on 
the other hand, supersymmetry is broken at tree 
level if and only if V has no zeroes. The predic- 
tions of strong coupling and weak coupling are 
quite different, and it is interesting to study the 
crossover from strong to weak coupling. 

3. MONTE CARLO SIMULATIONS 

We perform our simulations using the Green 
Function Monte Carlo (GFMC) algorithm [§). A 
discussion of GFMC in the context of the present 
problem can be found in Ref. Q; we only wish 
to remark the main features of the algorithm: 
the aim is to generate a stochastic representation 
of the ground-state wavcfunction, which is then 
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used to compute expectation values of observ- 
ables. Statistical fluctuations are reduced with 
the help of a guiding wavefunction, whose free pa- 
rameters are determined dynamically during the 
simulation. In order to keep the variance of ob- 
servables finite as the simulation proceeds, it is 
necessary to simulate a population of K walkers 
(field configurations at fixed time), and extrapo- 
late the results to if ^ oo. 

We focus on the case V — A20^ -I- \q4>; strong 
coupling always predicts supersymmetry break- 
ing; weak coupling predicts unbroken supersym- 
metry for Ao < 0; according to Ref. unbro- 
ken supersymmetry should be accompanied by a 
nonzero (0) (parity breaking). 

We presented preliminary results, obtained by 
GFMC, for intermediate couplings in Ref. our 
aim is to extend the study towards the continuum 
limit and to larger lattices. 

Perturbative computations show that 

A2 ~ aA^™, Ao ~ aAJi™ + aX'^""^ In (aM) , (1) 

27r 

where A^ is the adimensional lattice bare coupling, 
AJ''" is the renormalized (continuum) coupling, 
with dimension of m^, defined at the mass scale 
M, and a is the lattice spacing. We study, as 
A2 — > 0, the trajectory 

Ao = ^ ln(4A2), (2) 

ZTT 

corresponding to a perturbative RG trajectory 
(|l|); the effect of Ao is small in the range we con- 
sidered, therefore we expect Eq. (^ to be a rea- 
sonable approximation to a true RG trajectory. 

We estimate the correlation length from the ex- 
ponential decay of the connected correlation func- 
tion Gd ~ {4>n4'm)c averaged over all n,m pairs 
with |m — n| = d, excluding pairs for which m 
or n is closer to the border than (typically) 8. 
In our formulation, fermions are staggered and 
even/odd d correspond to different channels. 

We begin with the discussion of the case V = 
0.35 (/)^ -f 0.02, for which we have obtained the 
statistics of 4x10'^ GFMC iterations. The even-rf 
channel is plotted in Fig. Q it is very difficult to 
extract a correlation length, presumably because 
4> has a very small overlap with the lightest state 
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Figure 1. Connected correlation of (j) at even dis- 
tance for V = 0.35^^ -I- 0.02; the solid line is an 
exponential fit from distance 8 to 18; L = 34 and 
if = 100 data are consistent with the data shown. 



of the channel, and the value of ^ quoted in Fig. |l| 
should be considered tentative. The odd-c? chan- 
nel, plotted in Fig. ^, is much cleaner, and it is 
possible to estimate ^ with a good precision. 

For the other values of A2 , the situation is sim- 
ilar but with larger errors; we have a statistics of 
at least 10^ iterations, which we are increasing to 
4x10^. The values of ^odd follow nicely the ex- 
pected behavior ^ oc 1/ A2 , as shown in Fig. ||: the 
entire range 0.088 < A2 < 0.35 seem to be in the 
scaling region, with A2 = 0.5 a borderline case. 
The values of ^even have very large errors, and it 
is hard to draw any conclusion from them. 

We measure the ground state energy Eq along 
the trajectory (|^); the measurements have a very 
small statistical error, ranging from 1% for A2 = 
0.044 (where Eq/L w 10"^) to 0.1% for A2 = 0.5. 
We extrapolate to L 00 and K ^ oo fitting 
Eq/L to the form 

Eo/L = Soil + cl/L + ck/K). (3) 
Eq/L is plotted in Fig. ^: it seems to behave 
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1/^ = 0.165(4) 
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Figure 2. Connected correlation of at odd dis- 
tance for V = 0.35 + 0.02; the solid line is 
an exponential fit from distance 3 to 15; K = 100 
data are consistent with the K = 200 data shown. 
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Figure 3. Correlation length in the odd distance 
channel vs. A2; the solid line is a scaling fit. 



oc A2^^, while nai've scaling would predict oc X^. 
The value of Eq / L (disregarding this puzzling ex- 
ponent) and the lack of any signal for a break- 
down of parity (like a double-peaked distribution 
of (p) strongly hint that the trajectory (H) belongs 
to the phase with broken supersymmetry and zero 
{(j)). We are repeating the computation for tra- 
jectories with smaller Aq. 
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Figure 4. Ground-state energy density vs. A2. 
The dashed line is the result of a fit; the solid 
line shows the nai've scaling behavior. 



